Exact Thermodynamics of Disordered Impurities in Quantum Spin Chains 



A. Kliimper" and A. A. Zvyagin"-'' 
""Institut fiir Theoretische Physik, Universitdt zu Koln, Ziilpicher Str. 77, D-50937 Koln, Germany 
^B. I. Verkin Institute for Low Temperature Physics and Engineering 
of the National Academy of Sciences of Ukraine, 47 Lenin Ave., Kharkov 310164, Ukraine 

Exact results for the thermodynamic properties of ensembles of magnetic impurities with ran- 
domly distributed host-impurity couplings in the quantum antiferromagnetic Heisenberg model are 
presented. Exact calculations are done for arbitrary values of temperature and external magnetic 
field. We have shown that for strong disorder the quenching of the impurity moments is absent. For 
weak disorder the screening persists, but with the critical non-Fermi-liquid behaviors of the mag- 
netic susceptibility and specific heat. A comparison with the disordered Kondo effect experiments 
in dirty metallic alloys is performed. 
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In the last few years the interest in non-Fcrmi-liquid 
(NFL) behavior of magnetic metals and metallic alloys 
has grown considerably. A large class of conducting non- 
magnetic materials does not behave as usual Fermi liq- 
uids (FL) at low temperatures. One of the best known 
examples of such a behavior is the Kondo effect for multi- 
{n) channel electron systems: For an impurity spin less 
than ^ a NFL critical behavior results 0. The critical 
behavior of a single magnetic impurity can also be con- 
nected with a quadrupolar Kondo effect or non-magnetic 
two- channel Kondo effect [Q. However, for most of the 
dirty metals and alloys, in which the NFL behavior was 
observed, see e.g. the recent review [| andRefs. |,|, the 
magnetic susceptibility (x) and low temperature specific 
heat (C) usually show logarithmic or very weak power 
law behavior with temperature (T). The resistivity lin- 
early decreases with temperature showing a large resid- 
ual resistivity. The last property together with the al- 
loy nature of the compounds suggests that the disor- 
der (a random distribution of localized /-electrons or a 
random coupling to the conducting electron host) may 
play the main role in the low temperature NFL char- 
acter of such systems. The idea of (non-screened) lo- 
cal moments existing in disordered metallic systems has 
been already pointed out recently . It was proposed 
that near the metal-insulator transitions (or for the suf- 
ficiently alloyed systems far from the quantum critical 
point) disordered correlated metals contain localized mo- 
ments. The random distribution of their magnetic char- 
acteristics (i.e. their Kondo temperatures, Tk, which are 
characteristic energy scales for the crossover between the 
screened, or strong-coupling regime, and the weak cou- 
pling behavior) may be connected either with the ran- 
domness of exchange couplings of itinerant electrons with 
the local moments or with the randomness of the 
densities of conduction electron states Q . In Ref. Q] the 
results of the measurements of the magnetic susceptibil- 
ity, NMR Knight shift and low temperature specific heat 
have been reported. To explain the features of the be- 
havior it was necessary to assume weak disorder, with a 
Gaussian distribution of the Kondo temperatures. How- 



ever the model, which was used for the explanation of 
the experiment, was oversimplified Q: The magnetiza- 
tion of a single magnetic moment was approximated by 
the Brillouin function B{ah/T + bTx), where h is the 
external magnetic field, and a, b are constants. It was 
mentioned in Q that the data for the specific heat and 
Knight shift did not agree with the ones predicted by this 
simple theory, especially for nonzero values of the mag- 
netic field. The disagreement could be caused either by 
an inadequate representation of the Kondo magnetization 
by the single-impurity theory, or by the simple replace- 
ment T ^ T + bTx in the Brillouin function, or because 
of the not perfectly symmetric Gaussian distribution of 
the impurity couplings. The inhomogeneous magnetic 
susceptibility was confirmed very recently in |^ by muon 
spin rotation experiments. The role of the long-range 
(RKKY) coupling between the local moments was taken 
into account recently in Q (Griffiths phase theory). The 
latter gives a qualitatively similar behavior as for models 
with non- interacting local moments 

It is known that the behavior of a single magnetic im- 
purity in a one-dimensional (ID) antiferromagnetic (AF) 
Heisenberg spin S ~ ^ chain as well as the behavior of a 
single Kondo impurity in a 3D free electron host are de- 
scribed by similar Bethe ansatz theories ||l], [lo| , |ri]] , e.g. the 
magnetization and the low-temperature magnetic specific 
heat of the impurity for both models coincide. The spin- 
^ Heisenberg model is the seminal model for correlated 
many-body systems. Most of its static properties are ex- 
actly known. The spin-i magnetic impurity manifests 
the total Kondo screening with FL-like low temperature 
behavior of the magnetic susceptibility and specific heat. 
In other words, the moment of the impurity is quenched 
by the host spins, like the one by the spin degrees of free- 
dom of conduction electrons for the Kondo effect. On 
the other hand, for the integrable lattice models one can 
incorporate a finite concentration of magnetic impuri- 
ties 12 1 without destroying the exact solvability (it is 

where a single 
Hence, for 



impossible in the free electron host 
magnetic impurity can only be embedded), 
the random distribution of magnetic impurities we can 
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suppose that low dimensionality is not principal for the 
Kondo screening. The absence of the magnetic ordering 
in the NFL Kondo systems j|] also confirms this assump- 
tion and leads to the goal of our present investigation: To 
find exactly the thermodynamics of the disordered en- 
semble of spin-i magnetic impurities in the Heisenberg 
chain (with various random distributions of the impurity- 
host couplings) for arbitrary values of the external mag- 
netic field and temperature. It is the first study in which 
the thermodynamic characteristics of a disordered inter- 
acting many-body system are calculated exactly without 
any approximations. In this Letter we show that: (i) 
For several kinds of strong disorder of the impurity-host 
couplings the (Kondo) screening is absent; (ii) for weaker 
disorder the quenching persists, but with a NFL behavior 
of the magnetic characteristics. We performed a compar- 
ison of our exact results with previous approximate ones 
and with experimental data of NFL alloys. 

We investigate the thermodynamics of the quantum 
spin-^ AF chain with spin-^ impurities. The Hamilto- 
nian of the system has the form H = J2j + Himp, 

where the host part is Hj^j^i = SjSj-i-i (the host ex- 
change constant J is equated to unity). The impurities' 
part of the Hamiltonian has the standard form for the ex- 
actly solvable lattice Hamiltonians. Suppose we have an 
impurity distribution in which impurities are not nearest 
neighbors, then for the impurity, say situated between 
sites m and m + 1 of the host we obtain 



imp — Jimp^^m.imp ~^ -^rmp,m+l -^m,m+l 



(1) 



where Ji, 



1) ^, and [., .] denotes the commuta- 



tor. The coupling of the impurity to the host (Jimv) is de- 
termined by the constant 6. It was shown in ||Tl| , [l3|] that 
precisely this constant determines the effective Kondo 
temperature of the impurity via Tk oc exp(— 7r|0|): For 
energies higher than this crossover Kondo scale one has 
the asymptotically free impurity spin i, while for the 
lower energies the impurity spin is screened, and the 
usual FL-like behavior persists. In other words, 9 mea- 
sures the shift of the Kondo resonance of the impurity 
level with the host spin excitations, similar to the stan- 
dard picture of the Kondo effect in the electron host. We 
can independently incorporate any number of such impu- 
rities into the host chain, each of them will be character- 
ized by its own 6j. Hence we obtain an ensemble of the 
spin-i impurities with their own Kondo temperatures. 
The lattice Hamiltonian Eq. (|^) has additional terms, 
which renormalizc the coupling between the neighboring 
sites of the host, and three- spin terms. However it was 
shown in that in the long-wave limit such a lattice 
form of the impurity Hamiltonian yields the well-known 
form of the contact impurity-host interaction similar to 
the one of the usual Kondo problem [Q. The contact 
impurity coupling in this (conformal) limit is also deter- 



mined by the same constant 9. 

We managed to map our quantum Hamiltonian at fi- 
nite temperature to a classical system in 2D by means of a 
Trotter-Suzuki decomposition ||l5| . The geometry of this 
classical system is a square lattice with width L (^length 
of the quantum chain) and height N (=Trotter num- 
ber). The interactions on the lattice are four-spin inter- 
actions around faces with coupling parameters depend- 
ing on {NT)~^ and the interaction parameter 9i where 
i is the number of the column to which the considered 
face of the lattice belongs to. Note that the interactions 
are homogeneous in each column, but vary from column 
to column. We study this system in the limit N, L oo 
using an approach which is based on a transfer matrix de- 
scribing transfer in horizontal direction. The correspond- 
ing column-to-column transfer matrices are referred to as 
quantum transfer matrices (QTM). See Fig. ^ for an il- 
lustration of the model. 
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FIG. 1. Illustration of the geometry underlying the classi- 
cal model with four-spin interaction around faces and alter- 
nating coupling parameters from column to column. 

In general all QTMs corresponding to the L many 
columns are different. However, all these operators can 
be proven to commute pairwise. Therefore, the free en- 
ergy per lattice site of our system can be calculated from 
just the largest eigenvalues of the quantum transfer ma- 
trices (corresponding to just one eigenstate). For a dis- 
cussion of the homogeneous case see (l|Jl| and refer- 
ences therein. For our system we find the following set of 
non-linear integral equations for "energy density" func- 
tions of spinons a(x), a(x), A = 1 + a and A — 1 + d {x 
is the spectral parameter): 



k{x ~ y) In A{y) — k{x — y — in -\- ie) In A(y)] dy 



Inafa;) 



2T Tcoshx ' 



(2) 



with kernel function fc(x) — J duje^^^^^^^^^'^ / coshnuj. 
The corresponding equation for a{x) is obtained from 
Eq. (^) by exchanging i —i, h —h and a, A ^ d, A. 
The free energy per site / is given by 

^ ( . T f HA{y)Aiy)]dy 
fix)^eo{x)--J (3) 
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where cq is the groundstate energy. The free energy of 
the total chain with impurities is F = f{^Oj), where 
the sum is taken over all the sites (for sites without im- 
purities we get /(O)). These equations are easily solved 
numerically for arbitrary magnetic field values and tem- 
peratures. The random distribution of the values 6j (or 
of the Kondo temperatures for the impurities) can be de- 
scribed by a distribution function P{Oj). The details of 
calculations as well as the generalization of our model for 
the magnetically anisotropic case (important for systems 
with strong spin-orbital couplings ||]) will be reported 
later. It is worthwhile to emphasize here the simplicity of 
the derived equations: For each impurity there is only one 
parameter, the shift of the spectral parameter in the for- 
mula for the free energy per site Eq. (||). Then the exact 
solvability of the problem for any number of impurities 
permits to introduce the distribution of these shifts (or 
the strengths of the impurity-host couplings, i.e. the lo- 
cal Kondo temperatures). One has only two (non- linear) 
integral equations, Eqs. to solve, and the answer can 
be obtained for arbitrary temperature (in principle down 
to T = 10~^'*J) and magnetic field ranges. 

In Fig. H the results for (a) the magnetic susceptibility, 
and (b) the lincar-tcmpcrature coefficient of the specific 
heat 7 = C/T at zero field for a homogeneous Heisenberg 
chain, a single Kondo impurity, a Gaussian distribution 
of the host-impurity couplings, a Lorentzian distribution, 
and a so-called logarithmically normal JlSf distribution, 
which is characteristic for strong disorder, e.g. close to a 
critical point |^) are plotted as functions of the temper- 
ature. One can observe a clear qualitative difference be- 
tween strong disorder (Lorcntz and log-normal) and the 
other curves. It is clear that the divergent value of the 
magnetic susceptibility at low T for the strong disorder 
of the impurity-host couplings is connected with the fact 
that for most of the impurities their Kondo temperatures 
are lower than the temperature of the system. Therefore 
these impurities give rise to a Curie-like behavior of the 
susceptibility. This divergence disappears upon applying 
a finite external field which restors most of the FL-like 
behavior. 

Our results for the low-temperature thermodynamics 
are close to the ones of the perturbative calculations of 
random AF spin-i chains [p^|j2C| ]: 7 and x have weak 
power-law or logarithmic singularities. Our low tempera- 
ture results confirm the very weak dependence of the crit- 
ical exponents on the temperature. For the log-normal 
distribution we find critical exponents of 0.134 and 0.131 
for x(2^) and 7(T). In the case of the Lorentzian distri- 
bution these exponents are 0.846 and 0.730. The rather 
large deviations of these exponents are due to strong log- 
arithmic corrections in x(T) at low temperatures. For 
an illustration of this effect see also Fig. ^ showing the 
Wilson ratios 7/x which are non-universal, i.e. with NFL 
behavior and show infinite slope at T = 0. 
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FIG. 2. The temperature dependence of (a) x and (b) 7 at 
zero magnetic field for: Tlie liomogeneous Heisenberg chain 
(solid line); a single impurity (thin line); Gaussian distribu- 
tion of the parameter 9 (dashed line); Lorentzian distribution 
(iong-dashed hne); a log-normaf distribution (dashed-dotted 
fine) . The insets show log-log plots of the data. 

We like to note that former calculations were valid at 
most for low temperatures, while our approach is ap- 
plicable for any temperature and magnetic field scales. 
Furthermore, it is also known that the approximate re- 
sults [ p"9|j20[ incorrectly give zero or infinite susceptibility 
at r — > irrespective of the distribution, while its true 
value is finite for e.g. the homogeneous chain and the 
single impurity Our scheme, of course, perfectly 

describes the correct behavior. We also performed a com- 
parison of our results with the data of Ref . Q , with qual- 
itatively similar results. Note that any distribution P(A) 
of impurity couplings A(oc 1/6 in our parametrization) 
with finite P(0) corresponds to a distribution in 6 with 
Lorentzian tail. Hence the agreement of our data with 
|^,Q: infinite x and 7 for T — > 0. More clearly the com- 
parision with experiments can be seen from Refs. pl| , p2[ . 
By changing the concentration of impurities one can go 
from weak to strong disorder with quantative agreement 
with wide and narrow Gaussian and log-normal distribu- 
tions [El|. Critical exponents very close to those of our 
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Lorentzian distribution were observed in 121 
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FIG. 3. The temperature dependence of the Wilson ra- 
tio [1] at zero magnetic field. Depiction of distributions by 
lines as in Fig. ^. Note the infinite slope at T = due to 
logarithmic corrections. 

Generalizing model (1) by keeping exact solvability, it 
is possible to include (random) short- and long-range an- 
tiferromagnetic interactions of the special forms between 
the impurities themselves, see e.g. p^ , p3{ . However, 
these interactions do not affect the behavior of our disor- 
dered correlated spin system qualitatively, compared to 
the case without direct interaction between impurities. 
Here, a probability distribution P{9) with asymptotics 
\6\~" for large 6 leads to a divergence x(T) oc y-i/" 
which is weaker than that observed above. It was sug- 
gested in H] that the inclusion of other kinds of impurity- 
impurity couplings (of RKKY form, which violate the 
exact integrabilty) also do not change the behavior qual- 
itatively. However ferromagnetic impurity-impurity cou- 
plings can change the situation drastically, e.g. providing 
infinite x(T 0) even for weak disorder. The general- 
ization of our results to correlated electron systems with 
random impurities will be reported elsewhere. 

To conclude, we have constructed exactly the ther- 
modynamics of the Heisenberg antiferromagnetic spin- 
i chain with embedded disordered impurities. The re- 
sults are of high (numerical) accuracy and valid for ar- 
bitrary ranges of magnetic field and temperature. For 
strong disorder of the impurity-host couplings the local 
moments are non-quenched. For weak disorder of the 
host-impurity couplings, on the other hand, spin exci- 
tations of the host screen the impurities, but with the 
non-Fermi-liquid behaviors of the thermodynamic char- 
acteristics. The comparison of our theory with the data 
of a perturbative analysis and with those of magnetic ex- 
periments on disordered non-Fermi-liquids in rare-earth 
alloys shows qualitative agreement. 
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